Spatial solitons rays in periodic optical lattices 



a^ ■ 
o ■ 

o : 

Oh. 



> 
m 

cn 

o 
o 



R. Khomeriki^''^, J. Leon^ 
( ^ ) Laboratoire de Physique Theorique et Astroparticules 
CNRS-IN2P3-UMR5207, Umversite Montpellier 2, 34095 Montpellier (France) 
('^) Physics Department, Tbilisi State University, 0128 Tbilisi (Georgia) 

The light ray of a spatial soliton in an optical film whose refractive index is smoothly modulated 
(wavelength much larger than the typical soliton width) in both spatial directions is shown to 
possess chaotic regimes for which the propagation is erratic. This is interpreted as a parametric 
driven pendulum, obtained by a new perturbative approach of the Maxwell equation. These findings 
are then demonstrated to compare well to the eikonal law of light ray propagation (nonlinearity 
compensates diffraction). 



The spatial optical solitons, that result from a balance 
of nonlinearity and diffraction, manifest as stable self- 
focused light rays [H, S S S]- In a Kerr medium for 
instance, the paradigm model obtained from Maxwell's 
equations is the nonlinear Schrodinger (NLS) equation 
where the spatial direction of propagation plays a role of 
"time". Thus many standard classical and quantum tem- 
poral effects (as interaction, switching, tunnelling, etc.) 
can be usefully experimented on permanent regimes in 
the spatial domain 0, @, 0, S] ■ 

Recent works consider various spatial localization sce- 
narios in one and two dimensional optical lattices 
[lH , where the soliton localization width is comparable to 
(or even larger than) the intrinsic period of the grating. 
Much less studies have been devoted to the spatial soliton 
dynamics in smoothly inhomogeneous media, although a 
transversally modulated refractive index has been stud- 
ied in [3| ■ Our purpose here is to present the first 
discovery of erratic light rays that occur in a smoothly 
periodic optical lattice obtained as in [l4| by modula- 
tion of the refractive index in the two spatial directions. 
Moreover the modulation along the direction of propaga- 
tion will be assumed smoother than the transverse one 
leading then to a quasi two-dimensional problem. 

Then we propose a new approach to the spatial soli- 
ton description where the effect of lattice grating is con- 
sidered as a two-wave nonresonant interaction. We show 
that the chaotic motion is simply related to the paramet- 
ric driven pendulum which possesses periodic and chaotic 
regimes, instances of which are presented in fig[l]and ob- 
tained by changing the longitudinal modulation length 
of the optical lattice. Last but not least, this behavior 
is demonstrated to follow the rules of light ray optics 
(eikonal) which actually applies thanks to compensation 
of diffraction by nonlinearity. 

The Maxwell's equations in a non-magnetic film be- 
longing to the (x, 2;)-plane is written hereafter for the 
out of plane linearly polarized electric field component 
E{x, z, t) by assuming instantaneous response of the Kerr 
medium of susceptibility 
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FIG. 1: Graph (a): contour plots of the electric field envelope 
£{x,z) solution of NLS ^ in the potential Q with initial 
conditions ([5)l for two slightly different values of the injec- 
tion point xo = 460 for the red (r) ray and xq = 461 for the 
blue (b) one. The refractive index landscape Q is plotted for 
reference. Graph (b) is the driven pendulum approximation 
(|13p and graph (c)is the Eikonal model (|15|) . both simulated 
with the same initial conditions {A — 2 in all of those three 
cases). The lower graphs (d), (e) and (f) describe correspond- 
ing simulations with same initial conditions and parameters 
except for the longitudinal modulation Q replaced with 
(then here A = 0.5). 



The linear refractive index variations in the plane result 
from the optical lattice, thus n — n{x,z). The z-axis 
is chosen along the incident propagation direction and 
the standard approach to spatial solitons then consists 
in assuming slowly varying envelope, namely 



E = y^£{sx, £2^)e*['^*-'^'^l + c.c. + 0{e^ 



(2) 



g^ + Q;2)E~:^j^{n'E + xE'}^0. (1) 
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with wave number k = rigUj/c, where rig is a reference 
refractive index. The small expansion parameter e fol- 
lows from the smallness of the nonlinear Kerr coefficient. 
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thus one may simply define x = Moreover we are 
interested here in the physical situation when the spatial 
variations of the refractive index are also weak, namely 
rig — n{x, z) = C(£^)- In such a case the Maxwell equa- 
tion ^ in the leading order reduces to the NLS equation 
(e is scaled off) 



.d£ Id^e 



dz 



2dx^ 



n{x, z) 



f = 0, 



(3) 



which is written in dimensionless space variables by the 
use of new coordinates z ^ kz and x kx. By as- 
suming different scalings along the propagation direction 
z and transverse direction x, the grating index actually 
defines a quasi two-dimensional lattice, leading then to 
the above spatial f +1 dimensional equation in the exter- 
nal potential f — n/n^. 

So far nothing has been assumed for the relative di- 
mensions of the optical lattice wavelength v.s. spatial 
soliton extension. As a matter of fact preceding works 
as 0, [13, U I consider the case when the spatial soliton 
sees at least a few periods of the lattice while here we are 
interested in the case when the lattice modulation varies 
slowly with respect to the typical soliton size. 

To set an example let us assume a lattice periodic in 
both directions 

n{x, z) = n„ [l - xK s\n^{Kx) sin^(fiz)] . (4) 

with < xK ^ 1 such that — n{x, z) is always posi- 
tive. Then we choose initial conditions for ([3]) represent- 
ing the injection of a laser beam in the z-direction. 



E{x, 0) = asech [a^/x [x — xq)] 



(5) 



obtained from the exact soliton solution of ([3]) for n = n^, . 
The figdja) now displays two numerical simulations of 
© with 

a = 0.9, x = 0.01, K=0.2, K^^TT/2b, VL^Ky 

(6) 

and two very close initial injection points xo = 460 (red 
ray) and xq = 461 (blue ray). The two spatial soliton 
rays are plotted against the potential 1 — n/n„ for ref- 
erence. Therefore two very close initial positions result 
in very different trajectories indicating a chaotic nature 
of the soliton rays (note the trajectories destabilize al- 
though xq is in a region where the external potential 
vanishes). If now we increase VL, keeping all other pa- 
rameters values, we obtain a regime of a periodic stable 
trajectory presented in figlHd). 

We propose hereafter a comprehensive interpretation 
of such a behavior in terms of the equation of the para- 
metric driven pendulum. Our approach is to consider the 
Maxwell equation ^ where the light beam interacts with 
the index variations as a fundamental long-wave short- 
wave nonresonant interaction process. To that end we 



reconsider the reductive expansion method by following 
[lit and seek a solution as 

E = / J^*(^i, C)e*^(«^''^)e^['^*-'=^1 + c.c. + ■■■ , 
V 3x 

{n{x, z) - n„)/n„ = e'^V{£_2,0 

a =e[2:-</)(6,C)], ^2=e^x, C^e^z. (7) 

The property that the index variations are smooth with 
respect to the spatial soliton size is contained in the 
specific dependences on the slow variable ^ and very 
slow variable ^2- The resulting equation at order ap- 
pears with terms that depend solely on a single variable, 
^'(a, C) on the one side, 9{£,2, C) ^^id V{^2, C) on the other 
side. These terms thus decouple to eventually give 

.a* 152^- I ,2, „ de 

-'ac+w + '*' * = 5c = -^- 

It appears therefore that the effect of the optical lattice, 
through V{^2,C)> is not an external potential as in ([3|) 
but is contained in the coordinate drift (/)(a , C) defining 
£,1. Considering now (P) at order we get 



dC 



de 

'W2' 



(9) 



We eliminate first the phase e between (jH)) and ([9]), and 
second we turn back to initial spatial variables x and z 
to get 



d^cj) 1 dn{x,z) 
dz"^ n„ dx 



(10) 



Then the nature of equation ([5]) allows us to start with an 
explicit spatial soliton solution (in physical dimensions) 



ae 2'^ sech [a^(a; - (/))] 



(11) 



who se position Xs{z) is thus given by x^ — 4'{xg,z) = 
XjKibst. To compute the "acceleration" d^Xg/dz'^ we must 
remember that (j) is slowly varying in x and z and thus 
keep only the dominant orders, thus 



dz^ 



d'^4'{xs,z) 1 dn{xs,z) 



dz"^ 



dx.. 



(12) 



where the last relation follows from pO|) . This is the New- 
ton law in the "time" z for a particle of unit mass and po- 
sition Xs in the iime-dependent potential —n{xs,z)/ng. 

Considering now the example ^ used for the simula- 
tion in figUl the trajectory (|12p becomes the following 
parametrically driven pendulum model 1^ 17 1 



d^q 



^ -Asiii^t) sm{q), A = 2xK 



(13) 



where we have introduced a time t = Slz and angle 
q{t) = 2Kxs{z). These reduced units show that the main 
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control parameter is the above defined driving amplitude 
A. The parametric instability of the pendulum with ini- 
tial angle g(0) = and velocity q{0) = occur in the 
range 4/3 < A < 4, while the instability fastest grow- 
ing rate is reached for A = 16[(2/V3) — !]• However 
for different initial conditions, chaotic behavior may ap- 
pear at smaller (or larger) values of A [lil . IT7| . The 
graphs (b) and (e) in fig[T]show the result of simulations 
of (fT3|) with A = 2 and A — 0.5, respecively for initial 
velocity q{0) — and initial angles corresponding to the 
rescaled positions of the injected beams. The two differ- 
ent regimes presented result simply from a modification 
of the modulation period of the lattice along the propa- 
gation direction z. We observe indeed qualitative agree- 
ment, although trajectories should not be quantitatively 
compared since the regime is chaotic. 

In such solitonic regimes, nonlinearity compensates 
diffraction and the resulting light ray is a good candidate 
to geometric optics. As a matter of fact we demonstrate 
that the obtained soliton light ray equation can also 
be viewed as the eikonal of a light ray in a medium with 
variable index n(x, z). The eikonal equation written for 
the position f= (.t, z) in the plane reads 



d_ 

ds 



dr 
ds 



Vn, 



(14) 



where s is the arc length, namely ds^ — dx^ -\- dz^. We 
write this equation for the ray x = x{z) as we did for the 
soliton motion, and obain after some algebra that the 
above two components reduce to {x' = dx/dz) 



l\2 



x" = [1 + {x') 



1 



dn 
dx 



, dn 

dz 



(15) 



where the right hand side is eveluated on the trajectory 
X = x{z). To recover the soliton light ray, we must now 
use the assumption x' ^ 1, which is precisely the as- 
sumption made to obtain the NLS equation ([3|) from 
Maxwell equation llj. Taking also into account that 
n(x, z) 



n^ <^ ng one readily obtains from (|15p 



1 dn 

n„ dx ' 



(16) 



which is just the trajectory equation (|12p . Numer- 
ical simulations of have been performed for the 
same parameter values as precedingly (initial postion 
xq , vanishing initial veleocities and explicit expression of 
n{x{z),z)). The result is plotted in figmc) and figd^f), 
which confirms consistency with prediction of geometric 
optics. 

It is remarkable that such different models as NLS ([3]) 
in the external potential smoothly modulated in both 
spatial directions and the driven parametric pendulum 
concur to describe trajectory of a spatial soliton in 
a smoothly modulated (1-1-1) dimensional lattice. It is 



also striking that the soliton ray actually follows the rule 
of geometric optics in this context, which opens inter- 
esting perspectives on both theoretical and experimental 
aspects. 

One practical advantage of the parametric pendulum 
description is the prediction of the switching from regu- 
lar to chaotic (erratic) trajectories of the spatial soliton. 
This result might be useful for conceiving all optical ul- 
trasensitive noise amplifiers. Last, our result is expected 
to apply to the case of Bose-Einstein condensate with 
attractive nonlinearity in smooth space and time depen- 
dent optical lattices, for which a similar dynamics will 
occur but there in the time domain. 
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